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Clustering bifurcations are investigated by considering models of globally coupled map lattices. 
Typical classes of clustering bifurcations are revealed. The clustering bifurcation thresholds of 
the coupled system are closely related to the bifurcation structures of single map. In particular, 
cluster-doubling bifurcation induced period-doubling bifurcations and clustering induced chaos are 
found. At the onset of multiple-cluster states, equal-site-occupation-partition, and consequently, 
equal-phase-shift states [the so-called antiphase states reported previously, Pliys. Rev. Lett. 65, 
1749 (1990)] are always identified numerically. 


PACS numbers: 05.45.+b 

The investigation of globally coupled extended sys¬ 
tems has attracted a rapidly growing interest in recent 
years 1-6 . They arise naturally in studies of Joseph- 
son junctions arrays 1 , multimode laser 2 , charge-density 
wave 3 , oscillatory neuronal systems 4 and so on. A num¬ 
ber of intriguing and novel high-dimensional features 
have been revealed in these spatiotemporal systems. For 
instance, a curious and interesting dynamical state, so 
called antiphase state has been revealed both numeri¬ 
cally and experimentally 1,2 . Such a state is periodic in 
time, with each element of the system having precisely 
the same wave form. However, the motion of each ele¬ 
ment is just shifted by certain phase from its “neighbor”. 
This state is closely related to the clustering, which is nu¬ 
merically studied in the globally coupled chaotic maps 5 . 
The antiphase state is a clustering with equal occupation 
elements in each cluster. Up to date, the mechanism 
underlying this fascinating phenomenon has not been re¬ 
vealed. In this letter we will thoroughly analyze the bi¬ 
furcation mechanisms and phase diagrams of clusteriza¬ 
tions and reveal the interesting general features of equal- 
site-occupation-partition (ESOP) and equal-phase-shift 
(EPS) states at bifurcation points of clustering. (Here 
we call the antiphase state as EPS state because multi¬ 
ple phases may appear in clustering bifurcations.) We 
take the following globally coupled map lattice (GCML) 

L 

x n +i (i) = (1 - e)f[x n (i)) + - ^2 f[x n (j)\, i = 1, 

^ 3=1 

as our model, where n denotes the discrete time, i la¬ 
bels the lattice site with L system size. f(x) prescribes 
the local dynamics, and is chosen as the logistic map 
f(x) = ax(\ — x). e gives the long-range coupling 
strength. In Refs. 5, Kaneko et al presented very rich 
and interesting behaviors of (Q) for positive e. Negative 
e represents also many practical situations, such as an¬ 
tiferromagnetic coupling 6 and resistance coupling and so 
on. Therefore, it is useful to unify the investigations of 
Eq. ([[]) for both positive and negative e. 


First, we consider clustering bifurcations from the sim¬ 
plest spatially homogeneous configuration, so called co¬ 
herent state. After some simple algebra, the critical sta¬ 
bility condition of this coherent state can be explicitly 
shown as 

e c =l-e- Ao , (2) 

where Ao is the Lyapunov exponent of the single logis¬ 
tic map [ note, (12) is generally valid for any coherent 
state, whatever periodic or chaotic]. This critical stabil¬ 
ity boundary is shown in Fig. 1(a) with solid line. As 
e > e c , the coherent state always exists, and is locally sta¬ 
ble, while below the solid lines, the coherent state loses 
its stability and bifurcates to multi-cluster state. From 
(||) it is clear that, coherent periodic motions are always 
stable for positive e, they can lose stability only in the 
negative coupling regions. However, it should be empha¬ 
sized that many attractors may coexist with the coherent 
state in the regime above the solid line for large system 
size L and large a. 

A class of interesting states are multi-cluster states 
with ESOP (i.e., k clusters iVi = ••• = Nk, with 
Ni, i = 1, ■ • •, k being the occupation numbers of ith clus¬ 
ter), and in case of periodic motion, each cluster may 
have the same motion except some EPS. It will be shown 
that this kind of states appear naturally at bifurcation 
thresholds. Afterwards, a period-m state with k clusters 
(l)will be called TmCk state. It often happens that m = k, 
the evolving dynamics of the TkCk state can be much 
reduced as 

k 

x n+1 = {l-e)f(x n ) + ^-^2f(x j ). (3) 

3=1 

For k = 2, the solutions of Eq. (dj) read 


1 + a — ae ± W(1 + a — ae) 2 — 2(2 — e)(l + a — ae ) 

^ =- wT^T) -■ (4) 
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The stability conditions of the ESOP TkCk state can 
be analytically given by computing the products of Ja¬ 
cobi matrices of system (Jl]) . After certain matrix oper¬ 
ations, the stability analysis of system ( 0 ) can be much 
simplified to a block-diagonal form of the linear matrix 
of (y) as 


where 
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with = ax n +i{ 1 - X n+i ), i = 1,2, • • •, k, and x n+ k = 
k 

x n . M' = (1 — e) k a k x n (l — x n ) I, I is the (L — k) x 

n—0 

(L — k) unit matrix. Therefore, we obtain L — k de- 

k 


generate eigenvalues A = (1 — e) k a k x n (l — x n ), and 

n —0 

other k eigenvalues. If the absolute values of all eigenval¬ 
ues of J are less than one , the reference /c-cluster state 
is stable. For k = 2, The stability boundaries can be 
given explicitly. Increasing a from small value, the T2C2 
state appears from the spatially homogeneous state via 
saddle-node bifurcation (SN) and pitch-fork bifurcation 
(PK) for e > 0 and e < 0, respectively, at the following 
critical thresholds 


SN: a c 
PK : a c 


1 + 
2 + 
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1 — e ’ 


e > 0, 
e < 0, 
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and loses its stability via Hopf bifurcation (HF) and 
pitch-fork bifurcation at 


HF : a c = 
PK : a c = 


1 + 
1 + 


'1 + 


5 — 3e 


v (i~ e ) 2 : 
\J 2 + (e — 2) 2 


1 -e 


e > 0, 


e < 0, 


( 7 ) 


All these bifurcation lines and other clustering bifur¬ 
cation lines in the period-doubling region are shown in 
Fig. 1(b). It is remarkable that in Fig. 1(b) one can 
find a clear rule to describe the entire clustering bifurca¬ 
tion structure in e — a plane from the bifurcation points 
in the a axis at e = 0. Actually, all bifurcation points 
for a single cell (at e = 0) are multicodimension bifurca¬ 
tion points in the e — a parameter plane. For e > 0 one 
can find first order saddle-node bifurcation and second 
order Hopf bifurcation curves, while for e < 0 one can 


find second order pitch-fork bifurcation and first order 
subcritical bifurcation (SC) curves, all these bifurcation 
curves intersect with the a axis (e = 0) at the critical 
points for single cell. These beautiful bifurcation trees 
can be also found in chaos region associated with each 
periodic window. The enlarged regions of the rectangles 
in Fig. 1(a) for the bifurcations to ESOP three-cluster, 
and five-cluster states are shown in Fig.l (c) and (d), 
respectively. The bifurcation figures are similar to Fig. 
1(b). A difference of these figures from the two-cluster 
state is that for e < 0 these multi-cluster states appear 
from chaotic motions via saddle-node bifurcation rather 
than pitch-fork bifurcation. 

To give clear pictures about the clustering bifurcations 
we show the asymptotic states of system ([!]) in Figs. 
2(a) and (b), where all stable homogeneous states are 
plotted by diamonds, stable multicluster states by solid 
lines, unstable states (both homogeneous and inhomo¬ 
geneous) by dashed lines. The black regions represent 
stable quasiperiod motion. In Fig. 2(a) we fix e = 0.2, 
the system has only the coherent state below a c ~ 3.386 
(T1C1 for a < 3 and T2C1 for a > 3), the ESOP T2C2 
state occurs via a saddle-node bifurcation at a c , then the 
two states (T2C1 and T2C2) coexist in certain a interval. 
As a continuously increases, the coherent state undergoes 
a series of period doubling bifurcations leading to chaos, 
and then loses coherence at a ss 3.640. The T2C2 state 
subjects to Hopf bifurcation at a ~ 3.808. In Fig. 2(b) 
we fix e = —0.15, the bifurcations are essentially differ¬ 
ent from those of (a). The T1C1 state first undergoes a 
cluster-period-doubling bifurcation at a ~ 2.870 to cre¬ 
ate a stable T2C2 state. After a > 3.07, the coherent 
T2C1 state turns to be stable via subcritical bifurcation. 
In a large interval 3.4 > a > 3.07 the coherent T2C1 
state coexists with multicluster state, while the T2C2 
state undergoes further cluster-period-doubling bifurca¬ 
tion and Hopf bifurcation leading to chaos. At a ~ 3.4 
the T2C1 state bifurcates via cluster-period-doubling to 
form a T4C2 state. Fig. 2(b) is extremely interest¬ 
ing due to the following novel features. First, we find 
a cluster-doubling induced period-doubling. The value 
a ~ 2.870 is far below the period-doubling condition for 
a single map. Global coupling leads to cluster doubling 
at this parameter, that induces period doubling in time. 
Second, we find a cluster-doubling sequence 1-2-4 (and 
the induced period-doubling sequence). We expect that 
this cluster doubling cascade may proceed to a very large 
numbers of cluster. In our case this cascade is stopped 
at k = 4 by Hopf bifurcation at a ~ 3.316. Nevertheless, 
the tendency of cluster increasing bifurcations leading to 
chaos can be still seen in Fig. 2(c) for a < 3.4, where we 
plot the number of clusters vs. a for the states of Fig. 
2(b). It is found in Fig. 2(b) that chaos can appear for 
a < 3.4, where the nonlinear parameter a is far below 
the value for chaotic motion for the single map. In the 
same time in Fig. 2(c) we find the number of clusters di¬ 
verges (to the order of L ) in this chaos region. Then we 
conclude this chaos is made possible by clusterization, 
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and such is clustering induced chaos. It is remarkable 
that all analytical predictions in Fig. 1(b) are perfectly 
confirmed by numerical simulations of (□) in Figs. 2. 

In the above we focused on the discussion of ESOP and 
EPS multiple cluster states. On one hand, these states 
can be easily treated by analyzing Eqs. (§. On the 
other hand, these kinds of states appear generically and 
naturally under bifurcation conditions. For instance, in 
Fig. 2(a) around at a « 3.386 and from arbitrary initial 
conditions we can get T2C1 or T2C2 states. Whenever 
we get T2C2 state it must be an ESOP state (IVi = N 2 = 
^ if L is even, or Ni = -, N 2 = if L is odd). In 

Fig. 2(b) we run Eqs. ( 0 ) from random initial conditions 
at a = 2.8, then compute Eqs. (jlj) by gradually increasing 
a and by using the final state for the previous a as the 
initial state for the new a , we can surely get ESOP and 
ESP states for all cluster-doubling cascade. 

let us take a two-cluster state as an example. Usually, 
many two-cluster attractors with different Ni and N 2 
(TVi + N 2 = L) may coexist. However, at the onset of 
two-cluster state we always first find the ESOP state. In 
Fig. 3(a) and (b) we fix e = 0.2, L = 200, and plot the 
two-cluster states for the occupation numbers N\ = N 2 = 
100, and N\ = 98, N 2 = 102, respectively. It is really 
striking that in (b) by a slightly breaking the occupation 
balance the threshold for two-cluster state is considerably 
raised. In Fig. 3(c) we fix e = —0.3, L = 200, and plot 
N m vs. a, where N m is the possible smallest occupation 
number for the two-cluster states of Eqs. (1). Small 
noises are added to the system sites for wiping out glass 
states with extremely small basins. It is interesting that 
N m goes to at both left and right critical points. (This 
Phenomenon is also found by Kaneko et al for e > 0.) We 
find that the ESOP state naturally appears at both first 
and second order bifurcation points. 

Similar behavior happens also for general fc-cluster bi¬ 
furcations. In Figs. 4 and 5 we plot three and five cluster 
states, respectively, arising via saddle-node bifurcation 
from spatiotemporal chaos. In the two cases, the ESOP 
states arise much earlier than those with slight deviations 
from equal occupation partition for both e > 0 and e < 0. 
We have also examined many fc-cluster bifurcations and 
always find the same behavior. 

In conclusion we have revealed a detailed bifurcation 
structure for clusterization and found a cluster-doubling 
sequence and clustering induced spatiotemporal chaos. 
At any onset of fc-cluster state one always finds ESOP 
and EPS state. This is a general extension of the re¬ 
sults obtained by Wiesenfeld et al. However, in our case 
the natures of bifurcations and the links of clustering bi¬ 
furcations with the bifurcations of single site are clearly 
shown in Figs. 1 and 2, then these ESOP and EPS states 
of globally coupled systems can be predicted, based on 
the bifurcation structure of single site. 

In this letter we took the globally coupled map lattices 
as our models. However, the ideas can be extended to 
more general globally coupled systems. For instance, we 
have examined globally coupled Josephson Junctions, the 


clustering bifurcation features are qualitatively the same 
as those for the coupled map lattice systems. 
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FIG. 1 Bifurcation figure for homogeneous (coherent) 
state. Below solid curves the homogeneous state 
loses its stability, (b) Clustering bifurcations in 
period-doubling region. SN, HP, PK and SC rep¬ 
resent saddle-node, Hopf, pitch-fork and subcriti- 
cal bifurcations, respectively, (c) The blow up of 
the second rectangle of (a). The same as (b) with 
period-three window is considered, (d) The blow 
up of the first small rectangle of (a). The same as 
(c) with period-five window is considered. 

FIG. 2 Bifurcation sequences for certain couplings. Di¬ 
amonds, solid lines and dashed lines represent sta¬ 
ble coherent states, stable multiple cluster states 
and unstable states, respectively, (a) e = 0.2, at 
a « 3.386 a saddle-node bifurcation for two clus¬ 
ter state occurs, (b) e = —0.15, a cluster-doubling 
cascade induced period-doubling cascade leading to 
chaos is shown. Note the subcritical bifurcation for 
T2C1 state at a = 3.070, and the Hopf bifurca¬ 
tion at a = 3.316 which makes the cluster-doubling 
cascade incomplete, (c) The number of cluster, fc, 
plotted vs. a with e = —0.15. In the region of 
clustering induced chaos fc is of order L. 

FIG. 3 T2C2 states for L = 200. (a) e = 0.2, iVi = 
N 2 = 100. (b) e = 0.2, Ni = 98, N 2 = 102. T2C2 
state appears much earlier in (a) than in (b). (c) 
e = —0.3, L = 200. The number of the smallest 
site population N rn plotted vs. a. At the two bi¬ 
furcation points N m goes to -j, indicating ESOP 
state. 
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FIG. 4 Three-cluster states, L = 300. (a) e = 0.05, 
iVi, 2>3 = 100. (b) e = 0.05, Ah = 98, Ah = 

100, Ah = 102. (c) e = -0.1, Ah, 2,3 = 100. (d) 
(b) e = -0.1, Ni = 98, Ah = 100 ,N 3 = 102. At 
the onset of three-cluster state one can see only the 
ESOP state. 

FIG. 5 Five-cluster states, L = 500, the same behavior 
as in Fig. 4. (a) e = 0.005, Ah i2j 3 , 4,5 = 100. (b) 
e = 0.005, Ah = 94, Ah, 3 , 4 = 100 ,7V 5 = 106. (c) 
e = -0.005, All , 2 , 3 , 4 ,5 = i 00 . (d) e = -0.005, Ah = 
94, Ah, 3 ,4 = 100 ,Ah = 106. 
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